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ABSTRACT 



Context. Variability is a common characteristic of pre-main-sequence stars (PMS). Near-IR variability surveys of young stellar objects 
(YSOs) can probe stellar and circumstellar environments and provide information about the dynamics of the on going magnetic and 
accretion processes. Furthermore, variability can be used as a tool to uncover new cluster members in star formation regions. 
Aims. We hope to achieve the deepest near-IR variability study of YSOs targeting the p Ophiuchi cluster. 

Methods. Fourteen epochs of observations were obtained with the Wide Field Camera (WFCAM) at the UKIRT telescope scheduled 
in a manner that allowed the study of variability on timescales of days, months, and years. Statistical tools, such as the multi-band cross 
correlation index and the reduced chi-square, were used to disentangle signals of variability from noise. Variability characteristics are 
compared to existing models of YSOs in order to relate them to physical processes, and then used to select new candidate members 
of this star-forming region. 

Results. Variability in the near-IR is found to be present in 41% of the known population of p Ophiuchi recovered in our sample. 
The behaviours shown are several and can be associated with the existence of spots on the stellar surface, variations in circumstellar 
extinction, or changes in the geometry of an accretion disc. Using variability, a new population of objects has been uncovered that is 
believed to be part of the p Ophiuchi cluster. 

Key words. Stars: pre-main-sequence - stars: low-mass, brown dwarfs - stars: activity - stars: variables: general 



1. Introduction 



% 



Young stars have been known to be variable since iJovl d 1945b 
described the irregular behaviour of T Tauri stars. Photometric 
variability is thought to originate from several mechanisms re- 
lated to magn etic fields, ac cretion discs, and circumstellar ex- 
tinction (Herb st et al.lll994l) . Optical variability has been used 
as an excellent tool for characterising stellar and circumstel- 
lar environments of pre-main-sequence star s (e.gjGrankin et alj 
2007, 2 00J) and brow n dwarfs (e.g., ICaballero >et a?T 2004; 
Scholz & Eisloffel 2005). The mechanisms causing optical vari- 
ability in YSOs are also thought to be responsible for the occur- 
rence of photometric variations in the near-IR, especially suit- 
able for probing phenomena whi ch take place in the circum- 
stellar temperature regime (see Eiroa et al. 2002, and references 
therein). Furthermore, Kaas ( 1999) has shown that IR variability 
is a very useful tool for sorting the cluster population from the 
background field. For instance, it has the advantage of identify- 
ing young sources that do not show IR-excess or have small Ha 
emission, as is the case for weak-line T Tauri stars (WTTS), but 
are magnetically active dGrankin et al. 2008). The use of this rel- 
atively new technique should contribute significantly to obtain a 
full census of the young objects in star formation regions, one of 
the most important goals of galactic star formation studies. 

At a distance of 119±6 pc dLombardi et al.1 120081) . the 
p Ophiuchi cloud core is one of the nearest star forming re- 
gions, and therefore widely studied. Despite the proximity, few 
stars are optically observable due to the high visual extinction in 
the core, estimated to be 50-100 magnitudes ( Wilki ng & Ladal 



1983). That is also the case for most Myr old clusters, where 
much of the interesting population remains visibly obscured, 
and must be studied with IR observations. However, with many 
nearby star forming regions extending over degree scales and 
IR cameras usually having fields of only arc minutes, obtaining 
a single image of a region such as p Ophiuchi has been a ma- 
jor project until now, let alone attempting multi-epoch studies. 
This has changed with 2MASS which allows large area, though 
rather shallow, observations to be made, and near-IR var iabil- 
ity studies to become possible dCarpenter et alJ l2001. 2002). In 
a 2MASS JHK study centred near the Trapezium region of the 
Orion Nebula Cluster, where most star formatio n is thought to 
have occurred betw een 0.3 and 2 Myr ago dAli & Depoyll 1995c 
lHillenbrandlll997l) . ICarpenter et all d200ll) found that approxi- 
mately 45% of young stars were variable. 
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More recently, the WFCAM IR imager dCasali et al.ll2007l) 
has become operational on the UK Infrared Telescope (UKIRT). 
The large field of view (~0.8 deg 2 in four exposures) of this 
instrument has made deep variability studies possible for the 
first time in the IR. This paper presents the results of such a 
study, with multi-epoch H and K observations of the p Ophiuchi 
cluster used to search for variability on timescales of days, 
months and years. The study is enhanced b y including recentl y 
released IRAC/Spitzer data on p Ophiuchi dEvans et alj |2005). 
In Sect. 2, the observations and reductions for the large data 
sample are described. Sect. 3 explains the methods used in the 
search for variability together with the results of the analysis. 
These are discussed and explained in a physical context in 
Sect. 4, in an attempt to relate the variability behaviours ob- 
served to the possible physical causes. Conclusions and future 
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Table 1. Coordinates of the observations. 



Table 3. Photometric accuracy. 



Pointing RA 


Dec. 


1 16 27 20.0 

2 16 27 20.0 

3 16 28 18.2 

4 16 28 18.2 


-24 19 20 
-24 32 34 
-24 29 00 
-24 15 48 


Table 2. Log of the observations. 


UT date May 


June July 


2005 11,18,21,29 04. 


, 07, 14 


2006 24, 28^ 03. 


, 14, 23 04, 10 



Magnitude Photometric Accuracy 



H 



K 



Observations taken only for H band. 



prospects are given in Sect. 5. 



2. Observations 

The WFCAM on the UKIRT is a survey instrument developed 
at the UK Astronomical Technology Centre (UK ATC), pri- 
marily for the UKIRT Infrared Deep Sky Survey (UKIDSS) 
dLawrence et al . 2007). WFCAM is a wide field imaging cam- 
era operating in the near infrared fro m 0.83 pm to 2.37 u m in 
up to eight filters, including ZY JHK dHewettet al .1120061) . The 
camera uses four Rockwell Hawaii-II 2048 x 2048 18 /im-pixel 
array detectors with a pixel scale of 0.4". The four detectors are 
arranged in a 2 x 2 pattern and are separated by 94% of a de- 
tector width, for which four exposures are nee ded to survey a 
contiguous area of ~0.8 deg 2 dCasali et alj|2007l) . 

The data presented amount to a total of 14 epochs in the p 
Ophiuchi cluster. The observed region consists of four exposures 
which produce a final tile of ~0.8 deg 2 centred on the cloud's 
core. Each observing block was done in H and K bands where 
completeness limits of magnitudes 19.0 and 18.0, respectively, 
are achieved. TableQ]shows the central position (right ascension 
and declination) for each of the four exposures. The same region 
in the sky was observed randomly over scales of days, weeks, 
months and year sampling periods. Table |2]provides a list of the 
nights in which the observations occurred. 

To complement the study, Spitzer data from the c2d legacy 
project ( Evans et al.ll2003l) have been include d. The Ophiu chus 
molecular cloud has been mapped with IRAC dFazio et aU2 004) 
in the 3.6, 4.5, 5.8 and 8.0 /im bands, over a region of 8.0 deg 2 , 
which encompasses the WFCAM field. The data have been re- 
trieved from the prelimin ary c2d IRAC poin t-source catalogues 
of the third data delivery (Eva ns et al]|2005l) . 



2.1. Data processing 

WFCAM data are processed and archived by the VISTA 
Data Flow System Project, a collaboration between Queen 
Mary University of London, the Institute of Astronomy of 
the University of Cambridge and the Institute for Astronomy 
of the University of Edinburgh (Irwin et al., in preparation). 
The pipeline flat-fields the data, subtracts the counts from the 
background sky, detects and parameterises objects and per- 
forms the photometric and astrometric calibrations. The prod- 
ucts from the pipeline are a set of reduced uncalibrated indi- 
vidual exposures, and photometrically and astrometrically cal- 
ibrated stacked frames, as well as catalogues of sources de- 
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tected in the frames a nd ingested in t o the WFCAM Science 
Archive, described by lHamblv et alj d2008l) . A summary de- 
scription of the WFCAM pipeline and science archive is given 
bv lLawrence et al.1 d2007l) . 

Given the expected very low amplitudes of variation in 
young stars, a few ten ths or less of a magnitude in the near-IR 
(Carpente r et alj|2001l) . it is important to be aware of common 
problems intrinsic to the instrument and type of observations. 
To a large extent, many of these problems are eliminated by the 
pipeline during reduction and calibration. However, others do 
persist giving rise to false detections which show up in the cat- 
alogues. Although they represent a small fraction of the dataset, 
they need to be filtered out. A routine was written in IDL to im- 
plement an algorithm which finds and removes these spurious 
sources and implements the necessary constraints to ensure a re- 
liable final source list. These are briefly explained below. 

- Overlapping regions: the algorithm accounted for sources in 
the overlapping regions between detectors which have mul- 
tiple entries in the catalogue, keeping the source with the 
largest distance to the detector edge. 

- Magnitude limits: sets a cut-off limit at the brighter end 
(magnitudes 1 1.0 and 10.0 for H and K, respectively) due to 
saturation. At the fainter end, magnitudes 19.0 and 18.0 for 
H and K, respectively, mark the completeness limits which 
are above the 10cr level (see Table|3]for the average errors). 

- Object classification: objects flagged by the pipeline as stars 
were kept as well as objects which were classified as galax- 
ies but showed no difference in magnitude between different 
apertures. The motivation was to keep some nebulous YSOs 
while setting the criteria strictly enough to avoid galaxy con- 
tamination. 

- Cross-talk: electronic cross-talk occurs between the detector 
channels within a quadrant. It produces a sequen ce of spu- 
rious images in some or all the other 7 channels (Dy e et all 
2006). The algorithm removed sources located at up to seven 
multiples of 128 pix from bright stars. 

- Diffraction patterns: sources in the surroundings of saturated 
stars were removed since their photometry is not reliable. 
The radius of rejection was adjusted according to the bright- 
ness of the star. 

- Field edges: sources located on the outer field edge were re- 
moved since they were found to often have unreliable pho- 
tometry. 

Furthermore, the K band images from the observations made 
on 28 May 2006 were discarded since they were of lower quality 
due to technical problems. 

The catalogues for each night were merged for the two years 
separately. Sources were matched by requiring their position to 
be within an angular separation of 1 .0", since the position offsets 
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Fig. 1. Tile coverage in the p Ophiuchi cluster. The total area 
covered is 0.8 deg 2 . The spatial distribution of the detected 
sources reflects the cloud extinction. The visual extinction in 
the co re is estimated to be 50-100 magnitudes ( Wilkin g & Ladal 

1983). 



found range from 0.1" to 0.8". It was further required that each 
source had at least four epochs with detections in each filter. This 
criterion was implemented in order to reduce the probability of 
false transient detections and to reduce the effects of close dou- 
ble stars which are not resolved in all the nights, depending on 
seeing conditions. In total, 15 191 sources were detected across 
the two years of observations (Fig.|TJ. 

The internal accuracy of the photometry between nights was 
further improved by measuring and removing offsets, which 
can arise for several reasons. The photometric calibration of 
WFCAM data is done with 2MASS s tars in the observed frames 
dDveet all 120061: IWarren et alll2007l) . However, some of these 
stars could be variable, which would have an impact on the cali- 
bration. Furthermore, for regions in the sky with high extinction, 
as is the case for the central core of p Ophiuchi, the number of 
available 2MASS stars is reduced, which limits the photometric 
accuracy. Lastly, photometric offsets can also arise from small 
problems in flatfielding and sky subtraction from night to night. 
The removal of these offsets was done separately per detector 
and for each of the four sky pointings. Stars with magnitudes 
between 12.0 and 16.0 for K band and 13.0 to 17.0 in H band 
were used for the calculations, assuming that their magnitudes, 
in the median, do not change in time. The final offset for one de- 
tector at a given sky pointing was taken to be the average value 
of the differences in magnitude between the several nights and 
a reference night. This offset value, up to 0.03 magnitudes, was 
then removed from all the sources detected in that particular de- 
tector and sky pointing. 

Spitzer counterparts were obtained from the c2d IRAC point- 
source catalogues of the third data delivery, applying the rec- 
ommended criteria to select highly reliable samples, i.e. only 
considering sources which are not extended, as estimated by the 
source extraction, and with a quality flag of A or B in the four 
IRAC bands, which corres ponds to a signal- to-noise ratio ^10 
and k,7, respectively (see lEvans et al.ll2005l for a detailed ex- 
planation). Spitzer fluxes were converted to magnitudes using 



the IR AC zero magnitude flux densities defined by iReach et alj 
d2005l) . 



3. Results 

3. 1 . List of near-IR variable stars 

A variable star is one which displays photometric characteristics 
incompatible with photometric errors. The large sample being 
analysed required the use of statistical tools which can quan- 
titatively estimate the probability that the detected photometric 
variations are intrinsic to the star or its circumstellar material. 

The first tool used to detect variability was the reduced chi- 
squarex 2 (xl) of the magnitudes, computed using: 



1 



X v 



I 



(magi ~ mag) 2 



u- 



(1) 



where v is the number of degrees of freedom, N is the num- 
ber of measurements and cr, is the photometric uncertainty. In 
this case, the statistic represents the probability that the vari- 
ations result from Gaussian noise. Therefore, the reduced chi- 
square index entirely relies on the assumption that the noise is 
Gaussian, which is only approximately the case in any real set of 
observations in which there is variable seeing, bad pixels, detec- 
tor imperfections, etc. It is also very susceptible to outlier points 
and does not make use of correlated changes in multiband pho- 
tometry, further limiting its use. 

The second tool used takes advantage of the temporal co- 
herence of the star's light curve and also the correlation across 
H and K bands, allowing for true signals to be extracted from 
within the noise. The cross-correlation index (CCI) relies on 
simple correlation coefficients statistic s and was first pr oposed 
for use in the search of variable stars bv lWelch & Stetson] (11993). 
This method, which assumes that photometric errors for two 
different photometry lists should be uncorrelated, is effective 
as long as the difference in time between observations in each 
frame pair (at a given epoch) is small compared with the varia- 
tion period. The index is defined as: 



CCI 



N(N - 1) 



-5 



magy — mag 
cr H ,i 



t magf - 



mag 



(2) 



where N is the number of measurements and cr Hi , cr K j are 
the photometric uncertainties. 

A third index (CI//, CIk) was computed for each band sep- 
arately which takes advantage of the temporal coherence of the 
light curve and can be used to detect long term variations in a 
single band. It is defined as 



Ch 



1 



1 N(N - 1) 



s 



magf - mag\( magf +l - mag 



o~kj 



1 



0~K,i+l 



(3) 



for K band, and similarly for the H band. 

For each star, the;^ 2 ,, CCI, CI// and CI^-, were computed for 
each year separately. A star was classified as variable when it 
met one of the following criteria: 

- xl > 1^.3 in one or two bands, meaning only 1 to 2 (to 
be conservative) false detections in the total sample are ex- 
pected 
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Fig. 2. Cross correlation index (CCI). Variable stars show posi- 
tive values of the CCI greater than 2. 



- CCI, CIh or CIa: > 2.0. Numerical simulations have been 
used to study the behaviour of the correlation indices. For 
16000 stars (approximately the size of the dataset), the dis- 
tribution of values for the CCI due to random chance over 7 
epochs is symmetric about zero and has a maximum value 
of 1.8. The simulations for the individual band correlation 
indices showed the same behaviour, with a maximum value 
of the indices due to random chance of 1.9. Figure [2] shows 
the CCI as a function of H magnitude. Non-variable stars 
produce a probability distribution centred about zero due to 
noise while true variables show positive values of the CCI 
greater than 2 (to be conservative). All the correlation indices 
showed a similar behaviour. The width of the distribution 
of non-variable stars depends on the total number of epochs 
used in the calculation of the correlation indices, and so the 
threshold was scaled accordingly (CCI, CIh or CIa: > 2, 2.16, 
2.37, 2.65 for a star with detections in 7, 6, 5, or 4 epochs, 
respectively). 

The application of these criteria led to a raw list of 235 
candidates for variability. These were visually examined and 
59 objects were rejected, primarily because they turned out to 
be close double stars which fell inside the aperture, showing 
fluctuations in magnitude which originated from differences in 
seeing across the different epochs. Furthermore, faint objects 
(^T>15.5) which showed variability characteristics consistent 
with the properties of active field M dwarfs were also removed, 
since they are more likely to be affected by field contamination. 
Active low mass M -dwarf (dMe) stars hav e strong surface mag- 
netic fields dJohns-Krull & Valentil 1 1 996b and show significant 
levels of coronal activity with max imum amplitud e variations 
in the optical of <0.5 magnitudes (lBondar'1120021) . Assuming 
the same value of amplitude variation for IR variability (to 
be conservative), a simple starspot model dVrba et alJ [1986) 
predicts a maximum colour variation H-K of 0.015 magnitudes. 
All variables fainter than 15.5 in /T-band, magnitude amplitude 
variations less than 0.5 and colour variations less than 0.015, 
with a 3<t confidence, were therefore not considered in the 
final list of candidate members of p Ophiuchi since they were 
consistent with being active M dwarfs. This sample could, 



however, contain genuine cluster members, and so a table with 
the main properties has been included in an appendix. 

The final list contains 137 variable stars. Table [4] shows the 
photometric properties for each variable star, such as the coor- 
dinates, the average H and K magnitudes, the peak-to-peak am- 
plitude for H, K, and H-K, the slope in the variability colour- 
magnitude diagram (Sect. 3.4), the SED class as determined 
from IRAC/Spitzer colour-colour diagrams (Sect. 3.5), several 
variability flags (Sect. 3.4), and references for the known mem- 
bers of p Ophiuchi (Sect. 3.2). 

3.2. Known Population of the Ophiuchus Molecular Cloud 

Many imaging and spectroscopic studies have been done 
on the Ophiuchus molecular cloud across the spectrum, 
mainly at near- an d mid-infrared wavelengths, (e.g.,[V rba et al. 
1975; Elias 1978; Wilking etal. 1989; Comeronetal. 1993; 
Luhman & Rieke 1999; Bontemps et al. 2001; Allen et al. 2002; 



Natta et all 1 20061) but also in X-ravs dGagne et all 



Ozawa et al. 2005) and in the optical (e.g., iMartin et alJ 



2004 



1998 



Wilking et al. 2005). The WFCAM catalogue was merged with 
several catalogues of p Ophiuchi members (namely from the 



following studies. IComeron et al] 1993 ; [Bontemps et al. 2001; 
IGagne et all l2004t lOzawa et all I2005J Iwilking et al. 2005), 



which were chosen based on their accuracy to assess mem- 
bership, and also to better match the magnitude limits and 
spatial distribution of the WFCAM observations. WFCAM 
counter-parts are found for 128 previously known members of 
p Ophiuchi. Of the sources not merged with these catalogues, 
many are present in source tables of other surveys in the lit- 
erature, but they did not posses characteristics which enable 
them to be classified as young cloud members rather than back- 
ground stars or galaxies. A significant fraction of the p Ophiuchi 
population is not recovered simply because the objects are too 
bright, and therefore heavily saturated in the WFCAM images. 
Also the very extended objects are not present in the final cat- 
alogues since, to avoid galaxy contamination, this type of ob- 
ject was rejected. From the 128 known members with WFCAM 
counter-parts, 41% are variable stars, according to the previous 
definitions. In the sections that follow, variable stars which have 
been confirmed as members of the cloud by these studies are re- 
ferred to as members. Variables which have not been confirmed 
as members are referred to as candidate members. 



3.3. Magnitudes and colours of the variable stars 

The colours and magnitudes of the variable stars are an im- 
portant clue in investigating the youth and masses of these ob- 
jects. The colour-magnitude diagram, K versus H-K, shown in 
Fig. [3] displays the average magnitude and colour of the vari- 
able members of p Ophiuchi (open squares) and the candi- 
date members (filled circles). The age of obje cts in the core 
of p Ophiuchi has been found to be 0.3Myr dGreene & Meyerl 
1995; Luhman & Rieke 1999) and an extensive Ha survey by 
(Wilking et al. 2005) derived a median age for of 2. 1 Myr, with 
some members up to 3 Myr, for more widely distributed mem- 
bers. The curves in the diag ram show the model evolutionary 
tracks of the Lyon group dBaraflfe et alJ 119981 ; IChabrier et alJ 
2000; Ba raffeetalJl2003l) for a cluster at a distance of 119 
pc (Lombardiet al. 2008). The solid curve shows the theoreti- 
cal 1 Myr isochrone, for low mass stars ranging from 0.5 Mq 
down to 0.003 M0 (3My M „). The isochrone is a combination 
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Fig. 3. Colour-magnitude diagram for variable members (open 
squares) and variable candidate members (filled circles) of 
p Ophiuchi. Solid, dashed and dotted lines show the isochrones 
at 1 19 pc of the Lyon group (see text for detailed explanation). 
Lines parallel to the redening vector dRieke & Lebofskv| |l985) 
mark the separation between stars, brown dwarfs and planetary 
mass candidates for the different isochrones. 



oftheNextGenandt he DUSTY isochrones (iBaraffe et al.lll998t 
IChabrier et al.l l2000h. The dashed and dotted curves show the 
COND isochrones dBaraffe et all 120031) for 1 and 5 Myr, re- 
spectiv ely, down to 3My Mp . The lines parallel to the redening 
vector dRieke & Lebofskvll 19851) mark the separation between 
stars, brown dwarfs (M<0.075Mq) and planetary mass candi- 
dates (M<O.O12M0) for the different isochrones. The models 
indicate that the brighter candidates are above the hydrogen 
burning limit, but the faint sources would extend to very low 
masses. If an age of 1 Myr is considered, the faintest variable 
objects would have just a few Jupiter masses. However, when an 
older age is co nsidered (dotted line ), the objects would become 
more massive. Wil king et al.l d2005l) found a distributed popula- 
tion in Ophiuchus which is significantly older than that in the 
more highly extincted cloud core, with ages up to 3 Myr, so it 
is plausible that part of the variable objects are older than the 
core population. Although the CMD gives a good indication of 
the approximate masses of the variable objects, a spectroscopic 
follow-up is needed to investigate their nature. 

3.4. Characteristics of the variability 

The timescales of variability which can be studied with this 
dataset can be divided into two groups: intermediate and long 
time scales. The first applies to the stars which only show vari- 
ability in one of the two years of observations, and are likely to 
represent physical phenomena which have a maximum duration 
of a few months. The long term variables are stars that show 
variability behaviour across both years, which could be caused 
either by a single mechanism or be the result of repeated shorter 
term variability phenomena. In the list of member variables, 
77% of the stars show variability over the two years, while 
23% appear to be variable only in the first or second year. The 
values for the candidate members are slightly different, where 
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Fig. 4. AOC J162814.77-242322.5; variable candidate member 
star with no significant variations in the first year and large vari- 
ations in the second year. The panels show the H, K, H-K light 
curves. The left panels show the data points from 2005 May- 
June, and the right panels the photometry from 2006 May- July. 
The dashed lines show the mean values. 



52% of the objects are variable in both years and the remaining 
48% show variability in only one of the years of observations. 
Figure [4] is an example of a candidate member of p Ophiuchi, 
which shows no significant variations in the first year but large 
variations in the second year. Given the sparse cadence of the 
observations, a search for periods has not been attempted. 

The amplitude of the variations was measured as the 
peak-to-peak fluctuation in magnitudes and colour, for each 
year separately and also for the combined data. Figure [5] shows 
histograms of the peak-to-peak magnitude amplitude for all 
the variable stars, members and candidate members. They 
show a peaked distribution, with most of the amplitude values 
ranging from 0.01 to 0.8 magnitudes. Table [5] shows the mean 
and maximum values of the amplitudes for the members and 
candidate members. When the variation is measured across 
the two years, the average amplitude for the variable members 
is of the order of a few tenths of a magnitude. One extreme 
object, AOC J162636.8 1-241900.2, has been excluded from the 
calculations of Table [5] (see Sect. 4.2.5). The candidate member 
variables show mean K and H-K amplitudes of the same order 
of the variable members within l<x, and a mean H amplitude 
within 2<x. There is no statistical evidence that members and 
candidate members represent different populations. 

Most of the variable objects show a variation in H-K colour, 
which is on average 0.1 magnitudes but can be as large as 1 
magnitude (see Fig.|6]l. In some cases, this variation is correlated 
with the change in brightness, and can b e an important indicato r 
of the underlying causes of variability dCarpenter et al.l 12001). 
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Table 5. Amplitudes of the Variability 







Year 1 


Year 2 




Years 1 and 2 






Band 


Mean 


Max. 


Mean 


Max. 


Mean 


cr/ Vn - 1 


Max. 


Variable Candidates 


H 


0.19 


0.82 


0.20 


0.88 


0.27 


0.03 
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0.12 


0.70 


0.18 


0.02 
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Fig. 5. Histograms of the peak-to-peak magnitude amplitude of the variations in H, K and H-K. The top panels show the histograms 
for year 1 and the bottom panels for year 2. 



This behaviour is illustrated in Fig. [7] where the variations 
in the /T-band are linearly correlated with the changes in 
colour. The slopes in the colour-magnitude diagram K vs. H-K 
were computed using the linear regression method for objects 
which showed a highly significant linear correlation between 
change in colour and brightness, i.e. with < 1% probability that 
uncorrelated variables would yield a correlation coefficient at 
least as high as the assumed threshold. 36% of all the variables, 
members and candidate members, satisfy the linear correlation 
criterion and the respective slopes are listed in Table [4] The 
angle convention used by Carpenter et al. (2001) was adopted, 
where a slope angle of 0° represents a positive H-K colour 
change without a variation in K magnitude, and increases 
clockwise. Figure [8] is the histogram of the derived slopes 
which shows a binomial distribution with 16 objects having 
positive slopes, colour becoming redder as they fade, and 33 
showing negative slopes, colour becomes bluer as the star 



gets fainter. These distinct behaviours can be associated with 
physical processes and will be discussed later in this paper (see 
Sect. 4.2). 



3.5. Infrared excesses 

Young stars show infrared emission which originates from dusty 
envelopes and circumstellar discs surrounding the central object. 
lLada & Wilkingl £1984), based on the level of long wavelength 
excess with respect to a stellar photosphere emission, identi- 
fied three different classes which define an IR-excess or SED 
classification scheme: Class I, low mass protostars surrounded 
by an infalling envelope with large IR-excess; Class II, young 
stars with accretion discs and a moderate IR-excess (as classical 
T Tauri stars, CTTSs); and Class III, stars which no longer ac- 
crete matter from a circumstellar disc and show no IR-excess 
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Fig. 6. AOC J162814.73-242846.6; variable star with large vari- 
ations in magnitude as well as in colour. 
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(as weak-line T Tauri stars, WTTSs). The IRAC data from 
Spitzer allows the study of these objects in the mid-IR, where 
the excess contribution from discs and envelopes is predominant. 
The IRAC colour-colour diagram ([3.6]-[4.5] vs. [5.8]-[8.0]) 
was pr esented as a tool to separate young stars of different 
classes (lAllen e t al. 2004; Me geath et al.l2004l) . and was already 
used in the study of other star forming regions as, for exam- 
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Fig. 8. Histogram of the derived slopes in the K vs. H—K dia- 
gram. The slope can be used as a tool to distinguish the physical 
mechanisms behind the observed variations (Sect. 4.2). 



S ile, Taurus ( Hartmann et al.l2005t iLuhman et al.ll2006l). Se rpens 
Harvey et al.ll2007l) . or Chamaeleon II dPorras et alJl2007l) . 

The majority of populatio n of p Ophiu chi has been classified 
into SED classes. Bont emps et alj (1200 ll) used ISOCAM mid- 
IR bands (6.7 and 15.3 yum) to detect and classi fy 212 sources . 
Using ground-base mid-IR observations, Barsonv et al.l (120051) 
confirmed those results. Also in the X-ray regime, Chandra and 
XMM-Newton observations (Imanishi et al. 2001; Gagne et al. 
2004; Ozawa et al. 2005) have contributed to the classifica- 
tion. The variability catalogue contains 48 variable members of 
p Ophiuchi which were classified in these studies into Class I, 
II, or III. Figure|9]shows the IRAC colour-colour diagram for all 
variable objects with detections in the four bands, where the top 
diagram displays p Ophiuchi members and the bottom diagram 
the candidate members. The variable members are displayed ac- 
cording to the SED class assigned in the literature, i.e., Class I 
(open squares), Class II (cross es), and Class III (open triangles). 
The extinction vector is from Flahert vet al.l d2007l) . For com- 
parison, a sample of objects has been chosen from a region of 
the sky away from the central cloud with little extinction which 
should mainly contain field stars. 

The objects tend to cluster around three main regions of the 
diagram. Centred in the origin, [3.6]-[4.5],[5.8]-[8.0]=(0,0), 
are sources which have colours consistent with stellar photo- 
spheres and have no intrinsic IR-excess. These can be fore- 
ground and background stars, but also Class III stars which 
do not have significant circumstellar dust. In this region of the 
colour-colour plane, it is not possible to differentiate between 
young stars and contaminants. However, as it can be seen in the 
top panel of the diagram, many of these objects are confirmed 
young members of p Ophiuchi both from mid-IR and X-ray stud- 
ies. Another preferred region for obj ects in the diagram is lo- 
cated within the box defined bv lAllen et alj (120041) which repre- 
sents the colours expected from models of discs around young, 
lo w-mass stars. In fact, m any of the Class II objects classified 
bv lBontemps et al] (2001) lie within that range. However, some 
sources previously identified as Class II do not fall within the 
predicted limits. These objects have colours inconsistent with 
Class II sources (higher [3.6]-[4.5] colours than Class II but 
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lower [5.8]-[8.0] colours than Class I) and were previously clas- 
sified as candidate flat spectrum objects. Their location in this 
diagram does not confirm them as transition objects between 
Class I and II dBontemps et al.l200ll) since they can be explained 
as reddene d Class II sources. Finally, from models of infalling 
envelopes, lAllen et all d2004l) predicts the colours of Class I 
sources to have ([3.6]-[4.5]) > 0.8 and/or ([5.8]-[8.0]) > 1.1, 
which agr ees well with the two Class I objects identified with 
ISOCAM (Bontem ps et al.ll200ll) . Likewise, one of the sources 
classified with ISOCAM as Class III lies on the Class I/reddened 
Class II region of the colour-colour diagram. In the ISOCAM pa- 
per, this source was classified as Class III because it is located 
within the CS contours of the cloud, as defined by iLiseau et al.l 
d 19951) . According to the IRAC diagram, its previous classifica- 
tion is incorrect. 

From the 84 candidate members, 37 are detected in the four 
IRAC bands and plotted in the bottom panel of Figure [9] (filled 
circles). Four variable members are also detected in the four 



band but have no assigned SED class and therefore are also 
included (open diamonds). According to the criteria described 
above, they can be divided into Class I (2 objects), Class II (1 1 
objects), and possibly Class III (27 objects). 

The number of variables with detections in the 4.5 yum 
IRAC/Spitzer band is higher than for longer wavelengths, since 
this band has a higher sensitivity. Figure \W\ shows the H—K 
vs. AT-[4.5] colour-colour diagram for all variable objects with 
detections in the IRAC/Spitzer 4.5 //m band, where again the 
top diagram displays the members of p Ophiuchi and the bot- 
tom diagram the candidate members. The variable members 
are displayed according to their SED class from the literature 
(same conven t ion as in Fig. [9]). The reddening vector is from 
Flah erty et al.l d2007l) and the dash line follows the same redden- 
ing law. The top panel of Fig. [TUl shows a clean break between 
class I/II and class III stars (i.e., stars with and without circum- 
stellar dust emission) for ^T-[4.5]>1 and H-K>0.5. Applying 
this criterion to the variable candidate members (bottom panel), 
23 candidate members can be classified as ClassI/II, and 8 as 
possible class III. Furthermore, for objects present in both Fig. [9] 
and Fig. [TO] the SED classes determined are in agreement. The 
only exception are 3 variable candidate members which show 
a significant IR-excess in the /f-[4.5] colour but less at longer 
wavelengths. 

Combining WFCAM near-IR observations with 
IRAC/Spitzer mid-IR data, it is possible to classify 60% 
of the candidate members as Class I/II and possible Class III. 
Furthermore, the SED classes from the literature for the variable 
members are confirmed, with the exception of few objects (see 
discussion above). 



4. Discussion 

4.1. Membership 

The nature of the candidate member population is very impor- 
tant, since it determines the reliability of variability as an indi- 
cator of youth. In principle, some or all of the candidate mem- 
bers could be background objects - either variable stars within 
our own galaxy or extragalactic variable objects, principally 
Active Galactic Nuclei (AGN). An inspection of Fig. QT| how- 
ever, shows that this is unlikely to be the case. The candidate 
members do not appear to be distributed in the same way as the 
majority of background point sources against which the clouds 
can be seen - they are noticeably absent from the regions of low- 
est extinction for example, and appear to be more clustered. To 
put this on a quantitative basis th e 2-D Kolmogorov-Smirnov (K- 
S) test as outlined by Peacock] d 19831) was applied to candidate 
member and full point-source detection (non-variable) samples 
- the latter clearly dominated by the background population of 
stars and galaxies. The test is approximately distribution-free, 
which means the probability that the two samples are from the 
same population can be estimated. For 84 candidate member 
objects, a Z statisti c of 1.99 was de termined against the back- 
ground distribution. Peacock ( 1983J) gives a formula for calcu- 
lating the significance level, which shows that there is an 2% 
chance that the two are from the same population. This is an up- 
per limit since, for reasons already discussed, the variable can- 
didate members lie outside the regions of high extinction, as do 
the background sources. This introduces an unavoidable spatial 
bias, making the two distributions look similar. The K-S test sug- 
gests the two populations come from different distributions. Of 
course this is for the candidate population as a whole - there is 
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still the possibility that a subset are due to contamination from 
background variable objects. 

AGN are a potentially important source of contamination 
given the faint magnitude limits of the sample. The K - 18.0 
limit amounts to 0.08 mJy at 3.6 yum assuming an SED median 
index of -1.2 dBuchanan et a l. 2006). Treisteret al. (2006) use 
a Spitzer survey of the GOODS field to estimate AGN number 
density versus flux in the IRAC bands. Their data show that this 
flux limit should result in approximately 300 AGN detections 
per sq. degree. However, central to understanding contamina- 
tion by AGN i s the distribution a nd magnitude of extinction in 
p Ophiuchi. iKenvon et al.l ([1998), in a study of the reddening 
law, find that in the clear areas well outside the molecular gas 
contours in areas of least extinction, the H—K colour in most 
sources has only a small excess over intrinsic values - of or- 
der 0.2 in H—K. A similar excess is found in these data, corre- 
sponding to a visual extinction of 3-4 magnitudes. However, over 
much of the field the extinction is much higher. Therefore, for a 



proper evaluation, the number counts of K-ha.n& point sources 
(Fig. Q]) was used to calculate K band extinction in a grid of 
10 x 10 squares over the observed field. Th is was then used with 
the AGN cumulative flux distribution dTreister et alJ I2006T) to 
give the integrated number of AGN expected over the WFCAM 
field of view. After integrating over the extinction distribution 
a total of 33 AGN should be present in the p Ophiuchi field to 
the K band detection limit of 18.0. How many of these will be 
detected as variable? With the steep AGN cumulative flux distri- 
bution, more than 95% of the AGNs will lie in the faintest two 
magnitude bins, where the typical photometric error of the sam- 
ple is 0.05-0.1 magnitudes. Simulations with the 2-band cross- 
correlation index show that a clear variability detection with this 
index re quires an amplitud e of at least > 0.13 magnitudes in 
this case. lEnvaet alJ J2002I) show that the fraction of AGN with 
long-term IR variability larger than this is 28% or 9 objects in 
the p Ophiuchi field, compared to the candidate-member sample 
of 84. So while a precise number is difficult to calculate, it is 
likely that a small (but still non-negligible) fraction of the vari- 
able candidate-members can be explained as AGN. 

Another potential source of contamination is a popula- 
tion of background variable stars. The galactic model of 
Wain scoat et al.l (1 19921) was used as a useful tool in investigating 
this. Even though discrepancies between observations and this 
model exist, it is certainly accurate enough (factor of 2) for the 
purposes here. Schulth eis et al.l (|2000|) conducted an unbiased 
survey of a region in the galactic bulge, and found 720 sources 
variable in J and K band, representing 0.2% of the source de- 
tections in the region. These were mainly long-period variables 
- AGB stars at or beyond the tip of the RGB. A similar frac- 
tion of variables in the p Ophiuchi field (15191 sources) would 
predict 24 variables. However, the stellar populations are very 
different in the two fields. The bulge field is dominated by giants 
and includes many AGB stars; in fact the galactic model pre- 
dicts a number of AGB stars comparable to the total number of 
variables detected in the ISchultheis et alJ (12000) study. But the 
high galactic latitude of the p Ophiuchi field (+16.7) completely 
changes the result. Far fewer giants are expected, since the line 
of sight moves rapidly above the plane and bulge, so that the 
counts become dominated by the faint end of the dwarf lumi- 
nosity function. At this latitude, the model predicts < 1 AGB 
star/long-period variable in the whole field of view. It seems 
clear that the p Ophiuchi variable sample should be completely 
uncontaminated by these objects. 

It should be noted that in discussing the membership of any 
individual source, its location in the field is important, since 
any significant amount of extinction lowers the probability 
of background contamination substantially compared to the 
numbers calculated above. 



4.2. Understanding the Variability 

The characteristics of the variability observed in young stars in 
p Ophiuchi are very diverse, spanning large ranges in magnitude, 
colour and timescale. Furthermore, given the limited informa- 
tion and finite measurement errors, it is not possible to conclude 
with certainty which physical mechanism is causing variability 
in each star. In fact, in many cases the variability observed in a 
particular star could be caused by a variety of mechanisms. So 
the following approach was used. Known or plausible physical 
mechanisms were considered in turn as the cause of variability 
in each star. Mechanisms which were considered plausible ac- 
cording to certain criteria were then listed in Table |4] for each 
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star. In this way, an overall idea of whether the variability can be 
explained by conventional mechanisms, is obtained. A summary 
of the mechanisms, criteria and conclusions is as follows. 

4.2.1. Rotational modulation by cool starspots 

Cool starspots are known to exist both in CTTSs and WTTSs, 
and to cause variations in the stellar brightness. They are colder 
than the photosphere and arise from magnetic active re gions, 
analogously to solar sunspots dBouvier & B ertout 1989) . The 
fractional area coverage is of the order of 20% dBouvier et al.l 
1993; Her bst et all 1 19941) w hich determines the maximum am- 
plitud e of 0.4 magn i tudes dCarpenter et al.1 1200 U iHerbst et"aLl 
2002J). IHerbst et al.l d 19941) presented a classification scheme 
for optical variability, defining as type I variable stars with 
low, periodic, amplitude magnitude variations (a few tenths 
of a magnitude) caused by cool spots. Photometric variability 
studies of very low mass stars and brown dwarfs have identi- 
fied similar behaviours (e.g., Scholz & Eisloffel 2005). Using 
single-temperature blackbody models and input parameters in- 
ferred from optical st udies to model stellar spots behaviours, 
Carpenter et al.l d2001l) find that the expected amplitude colour 
variations produced by cool spots do not exceed 0.05 magni- 
tudes both in J—H and H—K colours. 

In the p Ophiuchi variable population, stars with amplitudes 
in H and K < 0.4 magnitudes and H—K colour variations < 0.05 
magnitudes show variations consistent with the existence of cool 
spots, as defined by the above parameters, and are classified as 
CS (cool spots; see Col. 7, Table|4]i. These amount to 17% of the 
variable members and 21% of the candidate members. Late-type 
field stars with magnetic activity could, however, still contami- 
nate this group of variables (Sect. 3.1). 

4.2.2. Rotational modulation by hot starspots 

Hot spots are interpreted as the impact points on the stellar 
surface from disc accretio n through magnetic field lines (e.g., 
ICalvet & H artmann 1992). They cover a smaller fraction of 
the stellar surface but the high temperatures can cause larger 



amplitude variations, which in the near-IR regime can be as 
high as 0.2-0.4 magnitudes for the JHK bands, and between 
0.05-0.12 in colours (Carpen ter et alj200l[). The times cales are 
shorter compared to cool spots (e.g jKenvon et al]|1994l) . In the 
IHerbst et al.l ([T994) classification scheme, a type II variability 
class is assigned to objects showing larger amplitude variations 
(which can be irregular or periodic), from short-lived hot spots. 
To ana lyse the behaviour produced by hot spots, Carp enter et al.l 
(2001)) use the information from the correlated colour and mag- 
nitude changes. The predicted slope from hot spot models has 
high, positive values, although given the simplicity of their mod- 
els, it is not possible to identify a specific range. 

Variable stars with amplitudes in H and K < 0.4 magnitudes, 
H and K colour variations 0.05-0.12 magnitudes, and a posi- 
tive slope in the colour-magnitude diagram (if a slope was de- 
termined) are classified as HS (hot spots). Hot spots can explain 
variability in 28% of the variable members and 30% of the can- 
didate members. 



4.2.3. Circumstellar Extinction 

Variations in extinction can originate from the intersection with 
the line-of-sight of infalling or orbiting material in the circum- 
stellar environment ( Skrutskie et al. 1996). Circumstellar extinc- 
tion variations will follow a reddening law if the grain size 
of the ma terial is comparable t o the interstellar grain size dis- 
tribution (Skrutskie et al.lll996l) . The interstellar redden i ng law 
E(H-K)/A K = 0.56 was adopted from lRieke & Lebofskvldl985l) . 
In the K vs. H—K diagram, this means that extinction variations 
can explain photometric slopes of 60°, with the object being 
redder when it is faint. The histogram of the colour-magnitude 
slopes (Fig. |6]l peaks at 45°-60°, showing that variations in ex- 
tinction can explain most of the observed positive slopes. 

Variable stars with positive slopes between 45°and 60°are 
classified as E (extinction). These represent 6% of the variable 
members and 4% of the variable candidate members. 



4.2.4. Accretion Discs 

Changes in the disc structure caused by mechanisms such as 
thermal instabilities, variable magnetic fields or warping insta- 
bility can lead to ne ar-IR variations on short timescales (see 
Carpenter et al. 2001, and references therein). Using simple disc 
models to study the consequences expected from geometric vari- 
ations in a circumstellar disc, Carpen ter et al.l (12001) find that 
the observed negative slope in the colour-magnitude correlation 
(Fig- ID can be explained by changes in t he inner disc struc ture. 
This scenario has also been suggested bv lEiroa et al.l (12002). 

Variable stars for which a negative colour-magnitude slope 
was found, where the star is bluer when faint, are classified as 
AD (accretion disc variations). In p Ophiuchi, 23% of the vari- 
able members and 25% of the variable candidate members have 
a negative slope. 

4.2.5. Extreme variability 

Although very rarely observed, another possible cause of 
variability rela ted to accretion disc s, is an Fu Orionis type 
outburst (e.g. IReipurth et al.l 120071 and references therein). 
AOC J162636.81-241900.2 (Fig. [IB has been previously clas- 
sified as a member of p Ophiuchi. Its light curve reveals a 
brightening of several magnitudes (3.5 magnitudes in the K- 
band) in approximately one year, in what could be a Fu Orionis 
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Fig. 12. AOC J162636.81-241900.2: variable member with a 
brightening of several magnitudes in both H and K bands. 



type burst. However, a minimum in H and K magnitudes of 
13.93 and 12.31, res pectively, have previously been reported by 
iBarsonv et al.l (1 19981) . meaning this object could instead be an 
EXor, which are stars that undergo a maximum in magnitude 
once every few years, thought to be caused by the massi ve in- 
fall o f circumstellar material onto the central star (see iHerbigl 
2008, and references therein). A spectroscopy study of this vari- 
able member (classified as EV) is needed to better understand its 
nature. 



4.2.6. Unexplained variability 

There are 32% of the members and 26% of the candidate mem- 
bers which do not fit into any of the criteria defined above. 
Among the members, 15 objects out of 17 have colours con- 
sistent with SED Class I/II stars. However the majority does 
not show a linear correlation in the colour-magnitude variabil- 
ity diagram. The most likely explanation is that several physical 
processes are taking place simultaneously, and the variability 
observed is the combination of all the intervenients. This sce- 
nario is further supported by the fact that 4 of these variable 
members show a positive slope which is either too shallow or 
to steep to be explained by changes in extinction, and colours 
inconsistent with hot spots rotational modulation. The same is 
true for the only candidate member with a determined, positive 
slope, which has also been classified as Class II, but does not 
fit into the variability categories. All the other candidate mem- 
bers with unexplained variability did not show a linear correla- 
tion in the colour-magnitude variability diagram, and very few 
have IRAC/Spitzer detections. The other characteristic which is 
valid for these objects, both members and candidate members, 
is that the mean H and K amplitude magnitudes, as well as the 
H—K colour variations, are twice as large as the ones seen among 
other variables. Although this is in conformity with the scenario 



just described, it also opens the possibility that more extreme 
variability behaviours are, after all, common among young stars. 
Only targeted photometric and spectroscopic monitoring of vari- 
able stars of different types, looking for a correlation between 
colours, magnitudes, and spectral features would help in isolat- 
ing the important physical mechanisms. 

5. Conclusions 

A multi-epoch, very deep near-IR survey was conducted in 
the Ophiuchus molecular cloud with the WFCAM/UKIRT 
Statistical methods, such as the reduced chi-square and correla- 
tion indices, were used in the search for variability, known to be a 
defining characteristic of young stellar objects. 137 variable ob- 
jects were found which show timescales of variation which can 
go from days to years, amplitude magnitude changes from a few 
tenths to 3 magnitudes, and significant colour variations. From 
the p Ophiuchi known population, 128 members have counter- 
parts in the WFCAM catalogue, and 41% are found to be vari- 
able. 

The dataset was merged with IRAC/Spitzer observations 
to further extend the information about the detected variables 
and using colour-colour diagrams, ([3.6]-[4.5] vs. [5.8]-[8.0]) 
and/or (H—K vs. A"- [4.5]), 50 variable candidate members were 
classified as candidate pre-main-sequence stars (23 Class I/II and 
27 possible Class III variables). 

The variability trends found were interpreted in a physical 
context by comparison with previous suggested classifications 
for variability in the optical and near-IR. The characteristics ob- 
served are consistent with the existence of cool or/and hot spots 
on the stellar surface, variations in circumstellar extinction, or 
structural variations in accretion discs. However, a large fraction 
of the variable population does not fit into the predicted param- 
eters for near-IR variability, which can be explained by the fact 
that the variability observed does not reflect a single event, but 
the net effect of several simultaneous processes. 

Finally, near-IR variability was used to discover a candi- 
date population of pre-main-sequence stars associated with the 
p Ophiuchi star forming region. The colour-magnitude diagram, 
K vs. H-K, provides evidence that many of the newly found low- 
luminosity variables have colours consistent with those expected 
from young, very low mass objects. To confirm the nature of 
these variable objects, a low resolution near-IR spectroscopic 
follow-up has been conducted for a small sample and results will 
be present in an accompanying paper. 
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AOC designation 1 


<H> Z 


<K> Z 


A// J 


A& 


AiH-K)' 


Slope 4 


SED 3 


Variability b/ 


Ref. s 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(degrees) 










J162555.67-240549.6 


13.97 


13.47 


0.05 


0.06 


0.04 






Y00Y 


CS 




J162558.29-241754.9 


17.23 


15.16 


0.45 


0.26 


0.50 






YYYY 






J162600.16-240427.3 


14.15 


13.71 


0.09 


0.09 


0.03 






Y0YY 


cs 




J 16260 1.37-242520.4 


12.30 


11.04 


0.45 


0.29 


0.17 


54 


II 


YYYY 


E 


2 


J162607.03-242724.1 


12.42 


10.40 


0.44 


0.31 


0.15 


65 


II 


YYYY 




2,3 


J162607.63-242741.4 


11.88 


10.27 


0.17 


0.16 


0.06 




III 


YYYY 


HS 


2,3 


J162609.92-243354.2 


11.77 


10.89 


0.21 


0.05 


0.17 




II 


YYY0 






J162615.01-244114.8 


15.19 


14.31 


0.31 


0.26 


0.22 






YYYY 






J162618. 13-243033.1 


12.98 


11.67 


0.06 


0.08 


0.04 


-53 


II 


YOYO 


CS,AD 




J162618.57-242951.4 


15.13 


13.55 


0.07 


0.08 


0.07 




II 


YOYO 


HS 




J162618.98-242414.0 


14.28 


12.15 


0.54 


0.47 


0.18 




II 


YYYY 




1,2 


J162619.10-240348.4 


11.42 


10.94 


0.11 


0.05 


0.12 




III 


Y00Y 


HS 




J162620.83-242839.5 


13.61 


11.87 


0.63 


0.47 


0.15 


59 


II 


YYYY 


E 




J162622.26-242407.0 


14.83 


13.69 


0.10 


0.14 


0.05 


-69 


II 


YYY0 


CS,AD 


1,2 


J 162623. 8 1-24 1829.0 


14.59 


13.49 


0.12 


0.04 


0.08 


20 


II 


YOYO 


HS 


1 


J 162624.29-24 1549.7 


15.81 


13.65 


0.14 


0.16 


0.03 




II 


YYYY 


CS 




J162624.41-244345.5 


16.44 


15.71 


0.23 


0.32 


0.37 


-36 




YYY0 


AD 




J162626.08-241 106.3 


15.20 


14.54 


0.20 


0.20 


0.11 






YYYY 


HS 




J162635.95-242058.7 


15.99 


13.59 


0.51 


0.53 


0.31 




II 


YYYY 






J162636.07-242404.2 


16.08 


13.18 


0.16 


0.13 


0.09 




II 


Y0YY 


HS 




J162636.81-241900.2 


14.86 


14.02 


2.78 


3.52 


1.16 


-72 


II 


YYYY 


AD,EV 


2,3 


J162637.41-240842.1 


11.21 


10.46 


0.14 


0.04 


0.13 




III 


YYYY 






J162637.78-242300.7 


12.60 


10.72 


0.77 


0.44 


0.34 




II 


YYYY 




2 


J162638.66-244144.4 


14.96 


14.15 


0.09 


0.10 


0.10 






Y00Y 


HS 




J162638.79-242322.8 


12.84 


11.52 


0.06 


0.07 


0.04 




II 


Y00Y 


CS 




J162639.92-242233.5 


14.69 


13.49 


0.08 


0.04 


0.06 




II 


YOOY 


HS 




J162642.14-243103.0 


13.68 


11.05 


0.88 


0.56 


0.36 




II 


YYYY 




2 


J 162642.5 1-24263 1.7 


15.63 


12.60 


0.06 


0.08 


0.05 




III 


YOYO 


HS 


1,2,3 


J162642.89-242259.1 


12.91 


11.37 


0.04 


0.06 


0.03 




II 


YOYO 


CS 


2,3 


J162643.90-240357.8 


16.28 


15.69 


0.19 


0.16 


0.33 


-24 




YOOY 


AD 




J162646.76-241908.4 


11.68 


10.29 


0.06 


0.07 


0.05 




III 


YOYO 


CS 


3 


J162647.05-244430.0 


11.20 


10.56 


0.10 


0.03 


0.11 




III 


YOYY 


HS 


2,4,5 


J162648.39-242834.8 


15.24 


12.46 


0.11 


0.09 


0.07 




II 


YOOY 


HS 


2,3 


J162648.48-242838.9 


14.07 


11.09 


0.49 


0.31 


0.21 




II 


YYYY 




1,2,3 


J162648.53-241227.5 


11.84 


10.55 


0.10 


0.08 


0.05 




III 


YOYY 


HS 




J162650.47-241352.3 


11.32 


10.81 


0.10 


0.04 


0.11 




III 


YOOY 


HS 


5 


J162651.52-241855.1 


17.39 


15.80 


0.52 


0.28 


0.47 




III 


YOYO 






J162651.97-243039.6 


16.55 


13.23 


0.08 


0.11 


0.06 


-57 


II 


YOYY 


AD 


1,2,3 


J162653.48-243236.2 


16.93 


13.15 


0.20 


0.20 


0.16 




III 


YYOY 




1,2,3 


J162654.77-242702.3 


14.97 


12.77 


0.21 


0.18 


0.11 




II 


YYYY 


HS 


2 


J162654.98-242229.7 


11.92 


10.11 


0.07 


0.06 


0.06 




III 


YOYY 


HS 


2,3 


J162656.27-241017.6 


14.91 


14.09 


0.12 


0.12 


0.04 






YOYY 


CS 




J162658.39-242130.2 


13.18 


11.37 


0.18 


0.20 


0.08 




II 


YYYY 


HS 


2,3 


J162658.64-241834.8 


13.13 


11.61 


0.05 


0.08 


0.04 


-63 


II 


YOYO 


CS,AD 


2,3 


J162658.65-242455.5 


17.02 


14.50 


0.22 


0.22 


0.08 




II 


YYYY 


HS 




J162659.04-243556.9 


13.87 


11.99 


0.38 


0.19 


0.20 


42 


II/III 


YYYY 




1,2 


J162702.99-242614.7 


15.83 


12.42 


0.25 


0.20 


0.11 




II 


YYYY 


HS 


2,3 


J162703.58-242005.5 


14.88 


13.51 


0.34 


0.29 


0.21 




II 


YYYY 




2 


J162704.10-242829.9 


13.12 


10.77 


0.19 


0.16 


0.15 




II 


YYYY 




1,2,3 


J162704.57-242715.6 


13.22 


11.10 


0.39 


0.28 


0.17 


47 


II 


YYYY 


E 


2,3 


J 162704.88-24442 1.2 


15.72 


14.90 


0.11 


0.08 


0.09 






YOYY 


HS 




J162705.66-244013.1 


16.81 


13.81 


0.14 


0.10 


0.09 




III 


YOYO 


HS 


2 


J162706.59-244148.9 


11.43 


10.71 


0.21 


0.11 


0.22 




II 


YYYY 




2,4,5 


J162706.76-243815.1 


14.19 


11.25 


0.38 


0.46 


0.13 


-70 


II 


YYYY 


AD 


1,2,4 


J162707.02-244338.9 


15.61 


14.61 


0.03 


0.04 


0.05 


-39 




OYOY 


AD 




J162709.34-244022.5 


13.72 


11.32 


0.08 


0.14 


0.07 


-61 


II 


YYYY 


AD 


2 


J162710.02-242913.3 


15.21 


14.16 


0.29 


0.21 


0.10 






YYYY 


HS 


1 


J16271 1.16-244046.7 


13.47 


10.55 


0.97 


0.62 


0.52 




II 


YYYY 




2 


J162711.71-243832.1 


15.10 


10.93 


0.77 


0.49 


0.36 


40 


II 


YYYY 




1,2,4 
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Table 4. continued. 



AOC designation 1 


<H> Z 


<K> Z 


A// J 


A& 


AiH-K)' 


Slope 4 


SED 3 


Variability b/ 


Ref. s 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(degrees) 










J162712.12-243449.1 


13.09 


11.32 


0.79 


0.55 


0.27 


58 


II 


YYYY 


E 


2,4 


J162714.34-243131.9 


16.65 


15.20 


0.41 


0.30 


0.15 


57 


II 


Y0YY 


E 




J162714.49-242646.0 


15.46 


12.08 


0.17 


0.25 


0.19 


-55 


I/II 


YYY0 


AD 




J162714.61-242312.6 


16.65 


15.14 


0.12 


0.22 


0.15 


-51 


III 


YYY0 


AD 




J162715.44-242639.8 


13.43 


10.55 


0.33 


0.25 


0.24 




II 


YYYY 




4 


J162715.50-243053.8 


16.49 


12.35 


0.19 


0.28 


0.16 


-53 


II/III 


YYYY 


AD 


2 


J162715.87-242514.1 


16.12 


13.15 


0.33 


0.27 


0.11 


61 


II 


YYYY 


HS 




J162717.19-243512.8 


16.59 


14.65 


0.14 


0.03 


0.15 




III 


Y00Y 






J162718.36-242426.3 


13.34 


11.20 


0.11 


0.08 


0.07 




III 


Y0YY 


HS 




J162718.37-243914.8 


15.60 


12.05 


0.17 


0.15 


0.07 




II 


YYYY 


HS 


1,2,4 


J162721. 78-242953.4 


15.18 


10.51 


0.64 


0.43 


0.41 


40 


I/II 


YYYY 




1,2,4 


J162721.81-244335.8 


13.40 


10.74 


0.06 


0.06 


0.03 




III 


Y0YY 


CS 


2,4 


J162724.63-242935.5 


14.99 


12.40 


0.12 


0.11 


0.05 




III 


YYYY 


CS 


1,4 


J162726.22-241923.1 


14.42 


12.91 


0.08 


0.09 


0.08 


-46 


II 


Y0YY 


AD 


2 


J 162726.57-242554.4 


12.33 


11.77 


0.18 


0.23 


0.10 


-66 


II 


YYYY 


AD 


5 


J162727.06-243217.6 


14.59 


12.30 


0.07 


0.07 


0.05 


-47 


III 


YOYO 


CS,AD 


4 


J162732.05-241956.3 


17.21 


15.47 


0.18 


0.06 


0.17 




III 


Y00Y 






J 162732.70-243242.4 


15.49 


13.00 


0.51 


0.37 


0.28 




II 


YYYY 




2 


J162732.84-243234.9 


12.86 


10.90 


0.12 


0.12 


0.09 




II 


YYYY 


HS 


2,4 


J162733.75-242234.9 


16.49 


15.66 


0.60 


0.38 


0.28 


52 


II 


YY0Y 


E 




J162735.58-240918.7 


14.85 


14.11 


0.05 


0.07 


0.04 


-62 




YOYO 


CS,AD 




J162736.14-240457.7 


12.97 


12.06 


0.28 


0.21 


0.09 




II 


YYYY 


HS 




J162737.23-244237.9 


14.77 


11.41 


0.18 


0.27 


0.17 


-58 


II 


YYYY 


AD 


2 


J162738.06-242527.8 


15.37 


12.70 


0.12 


0.04 


0.09 


20 


III 


YOOY 


HS 


2 


J162738.95-244020.7 


14.19 


12.42 


0.23 


0.35 


0.14 


-63 


II 


YYYY 


AD 


2,4 


J162739.88-243851.3 


15.41 


14.06 


0.05 


0.06 


0.04 




III 


YOYO 


CS 




J162740.28-240927.2 


14.47 


13.63 


0.04 


0.06 


0.04 


-56 


III 


YOYO 


CS,AD 




J 16274 1.60-244644. 8 


15.78 


13.84 


0.39 


0.15 


0.32 




II 


YYYY 




2 


J162741. 74-244336.2 


14.69 


11.88 


0.31 


0.22 


0.12 




II 


YYYY 




2 


J162744.98-240628.6 


11.78 


10.74 


0.66 


0.44 


0.47 




II/III 


YYYY 






J162747.08-244535.1 


12.94 


11.09 


0.08 


0.08 


0.04 




II 


YOYY 


CS 


2,4 


J162748.1 1-240155.9 


16.49 


15.90 


0.48 


0.38 


0.72 






OYYY 






J162748.24-244225.6 


16.92 


14.86 


0.54 


0.52 


0.18 




i/ii 


YYYY 






J162749.45-240422.2 


14.19 


13.61 


0.07 


0.05 


0.04 




in 


YOOY 


CS 




J16275 1.90-244629.7 


11.82 


10.45 


0.06 


0.06 


0.04 




in 


YOYY 


CS 


2,4 


J162800.20-240653.0 


13.64 


13.05 


0.06 


0.08 


0.03 






YOYY 


CS 




J162802.10-240844.9 


16.03 


15.38 


0.26 


0.25 


0.11 






YYYY 


HS 




J162804.54-243448.6 


15.11 


12.67 


0.11 


0.14 


0.08 


-54 


ii/iii 


YYOY 


AD 




J16281 1.08-242919.7 


11.26 


10.63 


0.18 


0.02 


0.18 




in 


YYYY 






J16281 1.46-241458.7 


15.88 


15.28 


0.10 


0.11 


0.07 






YOYO 


HS 




J162812.71-241 135.7 


12.09 


11.07 


0.14 


0.21 


0.14 


-45 


ii 


YOYY 


AD 




J162814.72-242846.5 


15.36 


14.99 


0.58 


0.65 


0.97 


-31 




YYYY 


AD 




J162814.76-242322.5 


11.18 


10.40 


0.34 


0.26 


0.09 




ii/iii 


YYYY 


HS 




J162816.41-240738.4 


14.24 


13.52 


0.13 


0.10 


0.10 






YYYY 


HS 




J162817.73-241705.1 


13.23 


12.30 


0.10 


0.15 


0.06 


-61 


ii 


YOYY 


AD 




J162819.82-240906.3 


12.42 


11.89 


0.30 


0.37 


0.47 


-35 


in 


YYYY 


AD 




J162821.73-241711.3 


17.20 


16.63 


0.28 


0.34 


0.47 


-33 




YOOY 


AD 




J 162822. 10-24 1229.6 


16.69 


15.96 


0.32 


0.33 


0.41 






OYYO 






J162824.38-242145.3 


14.93 


14.34 


0.11 


0.11 


0.03 




in 


YOYY 


CS 




J162838.34-242637.0 


11.12 


10.70 


0.13 


0.02 


0.14 




in 


YYOY 






J162841.53-243323.8 


11.13 


10.69 


0.08 


0.02 


0.09 




in 


YYOY 


HS 




J 162844.28-244724.7 


15.92 


15.42 


0.17 


0.16 


0.06 






YOYY 


HS 




J162845.88-240835.6 


13.98 


13.33 


0.22 


0.23 


0.15 




in 


YYYY 






J162847.02-242813.9 


12.74 


12.07 


0.63 


0.44 


0.21 


65 


ii 


YYYY 






J162847.38-242337.2 


14.17 


13.91 


0.07 


0.11 


0.05 


-65 




YOYO 


CS,AD 




J162847.61-242536.3 


12.92 


12.64 


0.07 


0.12 


0.05 


-64 


in 


YOYO 


AD 




J162848.69-243150.1 


14.34 


14.04 


0.23 


0.18 


0.09 






YYYY 


HS 




J162849.98-243518.1 


12.44 


10.85 


0.06 


0.05 


0.03 




in 


YOYY 


CS 




J162850.84-242832.3 


11.77 


11.46 


0.08 


0.02 


0.08 




in 


YOYO 


HS 




J162852.73-240926.2 


15.67 


15.01 


0.07 


0.09 


0.05 






YOYO 


HS 
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Table 4. continued. 



AOC designation 1 


<H> Z 


<K> Z 


AH* 


AK' 


A(H-tf) J 


Slope 4 


SED 3 


Variability 6 ' Ref. s 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(degrees) 








J162855.25-244515.6 


17.87 


17.33 


0.82 


0.36 


0.89 






YYYY 






J162856.53-244458.3 


17.75 


17.30 


0.53 


0.55 


0.72 


-35 




Y0YY 


AD 




J162856.79-243810.1 


11.16 


10.53 


0.13 


0.01 


0.14 


-4 


III 


YYYY 


AD 




J162856.94-243109.8 


11.51 


10.88 


0.07 


0.06 


0.04 




II 


YOYO 


CS 




J162857.78-242804.6 


11.08 


10.77 


0.12 


0.02 


0.13 




III 


YYYY 






J162904.32-242125.7 


16.49 


13.95 


0.10 


0.05 


0.06 




III 


Y00Y 


HS 




J16291 1.85-243955.9 


16.38 


15.32 


0.14 


0.10 


0.09 






Y0YY 


HS 




J162915.63-240958.5 


16.94 


16.26 


0.32 


0.25 


0.43 






YYYY 






J162916.32-243844.0 


14.04 


13.49 


0.30 


0.28 


0.03 


82 




YYYY 


CS 




J162919.97-240926.3 


14.69 


14.28 


0.05 


0.06 


0.02 


-67 




YOYO 


CS,AD 




J162923.28-240914.9 


13.97 


13.51 


0.05 


0.07 


0.07 




III 


YOYO 


HS 




J162924.58-240332.4 


14.48 


14.08 


0.14 


0.08 


0.07 


42 




Y00Y 


HS 




J162927.18-241020.0 


17.06 


16.61 


0.25 


0.45 


0.53 


-35 




YYOY 


AD 




J162930.31-244139.0 


11.25 


10.55 


0.08 


0.01 


0.09 


-7 


III 


YOOY 


AD 




J162931.01-243300.1 


15.47 


15.24 


0.12 


0.14 


0.15 






YOYO 






J162932.80-243305.9 


11.17 


10.90 


0.07 


0.06 


0.03 




III 


YOYY 


CS 




J162937.23-243426.8 


14.19 


13.90 


0.23 


0.25 


0.21 




III 


YYYY 






J162938.15-243618.4 


17.21 


16.87 


0.41 


0.32 


0.49 






YYOY 







J2000.0 IAU designation. 

Average magnitude over all available epochs. 

Peak-to-peak magnitude amplitude variation. 

Slopes in the colour-magnitude diagram A" vs. H-K. See Sect. 3.4 for details. 

SED class as defined from the combination of IRAC/Spitzer and WFCAM/UKIRT data. See Sect. 3.5 for details. 

First digit: star identified as variable from the CC indexes; second digit: star identified as variable from the xl statistics; third digit: star is 
variable in year 1 ; fourth digit: star is variable in year 2 (a star is a long term variable when both these digits are Y). When each of these conditions 
is verified it is denoted by Y (as in yes), otherwise the digits are set to 0. 

7 CP = cool spots; HT = hot spots; E = extinction; AC = accretion discs; FO = Fu Orionis. See Sect. 4.2 for details. 

8 p Ophiuch i mem ber s, according to the following studies: 1. IComeron et alj dl993h : 2. iBontemps et al.l d200lh : 3. iGagne et alj ( 120041) : 
4. lQzawa et all J2005h : 5. IWilking et alj i2005h . 
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Table A.l. List of variables with variability characteristics con- 
sistent with active field M dwarfs. The list could still contain 
young very low mass objects, possibly members of Ophiuchus. 



AOC designation 


<H> 


<A> 


AH 


AK 


A(H-K) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


J162554.64-244112.8 


18.20 


17.70 


0.51 


0.33 


0.52 


J162556.76-243001.1 


17.88 


16.30 


0.22 


0.23 


0.21 


J162609.95-240046.7 


17.22 


16.84 


0.10 


0.27 


0.32 


J16261 1.05-241208.8 


17.66 


17.18 


0.19 


0.40 


0.37 


J162612.86-241244.1 


17.71 


17.01 


0.22 


0.29 


0.16 


J162619.42-243233.8 


18.38 


16.98 


0.35 


0.31 


0.32 


J 16262 1.54-24 1242. 5 


18.61 


17.80 


0.60 


0.32 


0.48 


J162623.33-242351.4 


18.09 


16.47 


0.57 


0.30 


0.52 


J162631. 14-241208.5 


16.91 


15.65 


0.12 


0.11 


0.12 


J162647.88-240946.8 


16.42 


15.64 


0.27 


0.26 


0.16 


J162700.96-241230.6 


18.38 


17.72 


0.40 


0.30 


0.48 


J162705.45-240519.9 


17.02 


16.43 


0.20 


0.15 


0.13 


J162749.40-242919.3 


17.82 


16.61 


0.36 


0.20 


0.41 


J162749.80-243950.1 


16.94 


16.02 


0.39 


0.28 


0.18 


J162750.41-244042.7 


17.36 


16.43 


0.27 


0.07 


0.26 


J162800.15-241222.8 


18.35 


17.25 


0.58 


0.16 


0.53 


J162801.81-241305.0 


17.32 


16.44 


0.25 


0.12 


0.33 


J162803. 22-240536.7 


16.18 


15.63 


0.32 


0.20 


0.15 


J162814.41-244546.2 


16.88 


16.49 


0.21 


0.18 


0.08 


J162818.95-241319.1 


17.31 


16.68 


0.20 


0.26 


0.15 


J162840.66-243148.3 


17.72 


17.22 


0.42 


0.32 


0.35 


J 16284 1.89-242705. 6 


16.33 


16.02 


0.18 


0.20 


0.16 


J 162844. 20-240 115.1 


17.75 


17.21 


0.32 


0.28 


0.33 


J162846.90-240212.9 


16.82 


16.46 


0.08 


0.32 


0.28 


J162852.34-242849.1 


17.63 


17.37 


0.25 


0.25 


0.34 


J162858.30-243544.2 


17.50 


16.60 


0.22 


0.17 


0.08 


J162902.71-242837.7 


16.12 


15.78 


0.21 


0.13 


0.21 


J162905.65-242314.9 


17.79 


17.50 


0.22 


0.20 


0.36 


J16291 1.02-242240.0 


17.88 


17.50 


0.32 


0.20 


0.22 


J162919.04-241038.8 


17.72 


17.15 


0.27 


0.24 


0.22 


J162920.93-241936.8 


16.36 


15.90 


0.14 


0.11 


0.11 


J162929.44-243426.3 


18.00 


17.65 


0.57 


0.43 


0.59 


J162930.59-241907.8 


18.21 


17.68 


0.42 


0.43 


0.28 


J162933.26-243359.0 


17.96 


17.42 


0.47 


0.27 


0.59 


J162937.02-243035.2 


16.75 


16.40 


0.13 


0.19 


0.10 


J162939.40-242309.4 


17.19 


16.82 


0.30 


0.19 


0.21 


J162943.28-243056.4 


17.96 


17.57 


0.54 


0.52 


0.38 


J 162943. 49-243 102.7 


16.67 


16.45 


0.22 


0.18 


0.23 


J162945.00-243421.7 


17.71 


17.39 


0.46 


0.33 


0.28 



Appendix A: Other variable stars. 



